Lattice Calculations of 5 — )► K/K*l^l form factors 
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1 Introduction 



The exclusive rare semileptonic decays B — > Kl^l, B — )■ Kl^l, and the correspond- 
ing quark level process b — )■ sl'^l are mediated by flavor changing neutral currents 
(FCNC). In the standard model such processes arise at one-loop and are therefore 
I suppressed. Hence, these decays are good candidates to search for new physics beyond 

p ' the SM. Experimental measurements of -B — )• K/K*l^l~ decays have been performed 

D ■ by the BABAR, Belle, and CDF Collaborations P El E]- Most recently, new results 

^ have been reported by LHCb |11 El |6]. We expect to see increasingly accurate results 

from both LHCb and from planned high- intensity B factories 
^ \ To find evidence of new physics, it is necessary to compare the experimental results 

^ I with theoretical predictions from the SM. The current experimental results are con- 

^ ■ sistent with the SM predictions, but when the experimental error decreases, we need 

O ■ more accurate theoretical predictions. The accuracy of the theoretical prediction in 

B — )■ K/K*l^l~ is limited by the error of the hadronic matrix elements {B \ O \K/K*), 
which are parameterized by form factors. Some earlier theoretical calculations were 
based on the form factors calculated from Light Cone Sum Rules (LCSR). LCSR 
computes the form factors at low g^, where is the outgoing dilepton invariant mass 
squared. The form factors at large are extrapolated from low results, which ends 



. in a large error. Lattice QCD can calculate form factors directly at large q from first 



principles. Some earlier lattice QCD works calculated B — )■ K/K*l^l~ form fac- 
tors within quenched approximation P |9l HDl El [I2l [131 E]- These works obtained 
form factors at large q^ within 15%-20% accuracy. Moreover, the quenched calcula- 
tions have the systematic error from lack of sea quark, which is difficult to remove. 
Modern lattice QCD simulations include realistic sea quark effects, hence removing 
this systematic error. Recently, three lattice collaborations (FNAL/MILC, HPQCD, 
and Cambridge/W&M/Edinburgh) have started calculations of -B — )■ K/K*l^l~ 
form factors based on the 2+1 flavor ensembles generated by the MILC collabora- 
tion [T5]. The FNAL/MILC collaborations and HPQCD collaboration are working on 
the B Kl+l- decay HH US HE] . The Cambridge/W&M/Edinburgh collaboration 
are working on both of these two decays |T9l 120] . 
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Current methods in lattice QCD allow for calculations of 5 — > Kl^l^ with com- 
plete control over all sources of systematic error. This is not the case for B — )■ K*l^l~ , 
however, due to the fact that the K* meson is unstable, and needs to be treated as 
a resonance. Hence current lattice QCD calculations of this process contain an addi- 
tional systematic error that is difficult to quantify |19] . 

2 Lattice formalism 

The three new lattice QCD calculations use Nf=2+1 flavor gauge conflgurations gen- 
erated by the MILC Collaboration [21]. The MILC collaboration used the tree- level 
improved Liischer-Wise action for the gauge flelds and asqtad-improved staggered 
action for light sea quarks. These improvements suppress the lattice artifacts to the 
order of O^agO^) for the gluon fleld and 0{asa^), 0{a'^) for the fermion fleld \21\ . 
where a denotes the lattice spacing. These ensembles have four lattice spacings 
which are about 0.12fm, 0.09fm, 0.06fm, and 0.045fm. The FNAL/MILC collabo- 
rations employ ensembles on these four lattice spacings. The HPQCD collaboration 
and Cambridge/W&M/Edinburgh collaboration employ the ensembles with lattice 
spacings at a ~ 0.12fm and a ~ 0.09fm. Although these three groups use the simi- 
lar gauge ensembles, they employ different actions for valence b and s quarks. Both 
FNAL/MILC and Cambridge/W&M/Edinburgh group use the asqtad-improved stag- 
gered action for the valence strange quark. HPQCD uses HISQ action for valence 
strange quark [22]. The HISQ action is more improved than the asqtad action [22||2T]. 
For the heavy quark, the FNAL/MILC collaborations use the Sheikholeslami-Wohlert 
(SW) action [23] with the Fermilab interpretation for the b quark [21]. This action 
can be systematically improved and FNAL/MILC collaborations tune the b quark 
hopping parameter k and clover coefficient csw to remove the discretization errors 
through next-to-leading order {0{l/mb)) [25| [26] . The HPQCD collaboration and 
Cambridge/W&M/Edinburgh collaboration use (moving)-NRQCD method [271 [2H] 
for the b quark. The NRQCD method expands the relativistic QCD action by the 
order of Vb which is the velocity of the b quark. The heavy quark action is tuned 
to include the O^Aqfj-^/ml) corrections. Because these three groups are working on 
similar quantities with different discretization methods, it provides a good cross check 
for their results. 

3 The lattice QCD calculation on form factors 

The theoretical description of the B — )■ K/K*l^l^ process is based on the Operator 
Production Expansion (OPE). The low energy effective Hamiltonian for the b — )■ sl^l~ 
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transition is [301 [SH |32l [33] 



AC 

-Heff = -^\4,v;: ^c,(/i)o,(/i) (1) 

where OiS are four-fermion operators of dimension six. CjS are the corresponding 
Wilson coefficients. Most of the SM contribution is from the operators 07,9,10- The 
operator C7 is a photon dipole operator and C9,io are semileptonic operators. Theo- 
retical predictions are calculated from 'Heff and contain hadron matrix elements that 
are parametrized by form factors. For B — )■ Kl~^l~ decay, there are three form factors 
/+, /o, and /t 

{KUsrblB) = Uiq') (^p-,+p-^ - ~ ^^ g^) + /o(g^) ^^ ~ (2) 
{K\zsa^-^b\B) = -l^^^^p'^k'^-p-'knqu, (3) 

where pb and px are the B meson and kaon momenta, respectively. FNAL/MILC 
collaborations and HPQCD calculate these two matrix elements in the B meson rest 
frame. The q"^ becomes + — 2MbEk in this reference frame. The matrix 
elements are reparametrized as 



{K\i-s-f>^h\B) = ^/2M~B[v>'f\\{EK)+p^JdEK)\, (4) 

where = p'^/Mb is the four-velocity of the B meson and = — [px ■ v)v^. 
The form factors f\\ and are solved from the temporal and spatial components 
of the matrix element of the vector current. Finally, the form factors (/+, /o) are 
reconstructed from /y and by 

A = ^^?^J(MB-Ef:)!t + (Ei-Ml)f^]. (6) 

A similar method can be applied to fx'- 

^ MB + MK{K\tba''s\B) 

For B — )■ K*l^l^ decay, there are four hadronic matrix elements related to the 
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theoretical predictions. They are 

{K*{k,e)\sYb\B{p)) = M^y^/"'"'<KP'. (8) 
{K*{k,e)\s^^^^h\B{p)) = 2MK^Aoiq'f—^q'' + iMB + MK*)A,iq^)ie*''-'—^q'^) 

q" {K*{k,e)\sa^,b\B{p)) = m{q^)e^u.pe'pPnk. (10) 
q^ {K*{k,e)\-sa,,^,h\B{p)) = 2iUq^)[el{Ml - M],,) - {e* ■ q){p + k) ,] 

2iT,{q^){e* ■ q)[q, - J^T^^iv + k),] (11) 

Like the f+fl^T case, lattice QCD calculates the form factors by calculating a particular 
component of the matrix elements. Details on B ^ K*l^l^ form factors are in 

Ref. uniEni. 



4 Preliminary results on B — > K/K*l^l form fac- 
tors 

4.1 Extract form factors on the lattices 

First of all, we extract the meson masses and energies from the two-point correlation 
functions measured on the gauge configurations: 



c,{t-p-) = j2{op{x,t)ouo,o))e--^ = E(-i)"^^ ' 



:i2) 



2Ep 



where P represents the meson we want to study and Cp is the interpolating operator. 
If we insert a complete set of states, the two-point correlation functions are decom- 
posed into the contributions from different energy levels. The m labels the complete 
set of states that contribute to the sum. The factor (—1)™* arises only with staggered 
valence quarks. We are interested in the ground state energy and the excited states 
contributions can be safely neglected at large enough t. The FNAL/MILC collabora- 
tions use different smearing functions for B and K mesons to improve the accuracy 
of the results [SU [16]. The Cambridge/W&M/Edinburgh group applies the all-to- 
all propagator technique [35] to improve the signal-to-noise ratio. The ground state 
mass and energy in the current calculations are well- determined and have sub-percent 
statistical errors. 
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Lattice QCD extracts hadronic matrix elements from the three-points correlation 
function. For example, to determine /+ and fo in B Kl^l^ decay, we measure 
C3^fj^{t,T;'pK) which is defined as 

Cs,,it,T;pK) = 5^e*^--^^(O,,(0,0)V;(t,y)OUT,f)) (13) 

where V^=is'yfj_b. T denotes the location of the sink operator. Similar to the two- 
point correlation function, if we insert a complete set of states to the three-points 
correlation function Cs^^, it is decomposed into a sum of energy levels as 

C3,^(t,T;p) = ^(_i)-t(_i)-(r-*)A--e-<'*e-<'(^-*) , (14) 



where 



(15) 



The y4|; contains the matrix elements we want. The three-points correlation function 
also has the contributions from the excited states and opposite-parity states. Different 
methods are used to extract A^^. The FNAL/MILC collaborations apply an iterative 
averaging trick [34J to suppress the oscillating states contributions. The HPQCD 
and Cambridge/W&M/Edinburgh group fit the two-point and three-point correlation 
functions simultaneously to extract A^. They use the constrained fit technique [361 
, which helps to resolve the information from excited states. 



4.2 Chiral-continuum extrapolations and z-expansion fit 

Although the strange quark mass (m^) used in numerical simulations is typically 
close to its physical value, the light {u, d) quark masses in current simulations are 
usually larger than their physical values. Very recent new simulations include light 
quarks with masses at their physical values |3S1 SHI IH] , but they have not yet been 
used to analyze B meson decays. The FNAL/MILC collaborations perform combined 
continuum-chiral extrapolations using HMSxPT. 5*^7(3) HMSxPT was tested in B 
and D semileptonic decays |3H |12]. Some preliminary results from FNAL/MILC 
collaborations suggest the SU{2) HMS^PT might be a better effective theory in the 
B — ?► Kl^l^ process. The continuum form factors from lattice data are at the small 
Ek regime, because the discretization error becomes larger at large Ek in lattice 
calculations. Moreover, the correlation functions become noisier and HMS^PT is not 
reliable at large Ek- To have the form factors on the whole range, the FNAL/MILC 
collaborations use the z-expansion to extrapolate the lattice results to low range. 
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The z-expansion fit is a model-independent parametrization of form factors on the 
whole range. It maps the variable of to z: 



where t± = {Mb ± Mk/k*Y- to is selected to keep the absolute value of z smaller 
than 1. The form factors are then parametrized as: 



where B{q^) = z{q'^,mj^) is called the Blaschke factor, mji denotes the location 
of the pole in form factors, ipiq^ito) is called the outer function. In fits of lattice 
(or experimental) data, the 2;-expansion is truncated at some finite order. Different 
choices for B and yield different expansion coefficients. Generally, can be chosen 
to keep the small and hence the truncation error is well controlled. In addition, 
one can obtain bounds on the a, based on unitarity and heavy quark power counting. 
Hence the ^-expansion provides us with a systematically improvable description of 
the q^ dependence of the form factors. The preliminary results on form factors in 
B — > Kl^l~ from the FNAL/MILC collaborations are summarized in Fig. [H The 
continuum form factors at large (g^ > 15GeV^) are obtained from chiral-continuum 
extrapolations with SU{2) HMSxPT. The form factors at low q^ (g^ < 15GeV^) are 
from the 2;-expansion fit. The systematic errors from chiral-continuum extrapolations, 
heavy and light quark discretization, renormalization factors, scale determination, 
light quark mass determination, and finite volume effect are included. The total error 
at large g^ is about 5%, which is more accurate than the previous quenched results. 

The HPQCD collaboration and Cambridge/W&M/Edinburgh collaborations use 
modified 2;-expansion in the extrapolation of lattice form factors to continuum. As 
described in the last paragraph, the 2;-expansion fit provides a way to parametrize 
form factors on the whole g^ range. The HPQCD collaboration uses a modified version 
of z-expansion, where heuristic discretization and light quark mass dependent terms 
are added in order to perform combined chiral, continuum and shape fits [13]. The 
HPQCD collaboration applied it on the D semileptonic decays [13] and is planning 
to use this method for the B — )■ Kl^l^ process. The Cambridge/W&M/Edinburgh 
collaborations also use the same method for the extrapolations in S — > K/K*l^l^ 
decays. The preliminary results from these two groups are summarized in Fig. [2l 
The left panel is an example for /+,o,t in the B — )■ Kl^l~ process from the HPQCD 
collaboration. The HPQCD collaboration measures form factors with 1% accuracy at 
large g^ on the a ^ 0.12fm MILC lattice ensembles (only statistical error is included 
here.). More measurements will be done in the future on the a ^ 0.09fm MILC 
ensembles. Similarly, the right panel is the preliminary result of the Ti_2 in i? — t- 
K*l^l~ process from the Cambridge/W&M/Edinburgh collaborations. 




(16) 




(17) 
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vs. f+ and fg z-expansion fit 



vs. fj z-expansion fit 




Figure 1: Preliminary result of /+.o (left) and /r (right) on the whole range from 
FNAL/MILC collaborations. They obtain the form factors directly in the range > 
15GeV^ , and extrapolate them using the z-expansion into the region of < 15GeV^. 

5 Summary 

Lattice QCD can calculate the form factors m. B K/K*l^l~ decays from first prin- 
ciples. Previous quenched calculations obtained the form factors with an uncertainty 
of 15%-20% and over a range of q^ that was limited to large values. The FNAL/MILC 
collaborations, HPQCD collaboration, and Cambridge/W&M/Edinburgh collabora- 
tions are working on the new calculations with dynamical QCD configurations, which 
include the sea quark effect. The analysis techniques like constrained fit, SU{2) 
HMSxPT, and modified 2;-expansion are used in these new calculations. The pre- 
liminary B — )• Kl^l~ form factors from FNAL/MILC collaborations have a total 
(combined statistical plus systematic) accuracy of about 5% at large q^. Their final 
results will include a detailed and complete systematic error budget. We can expect 
increasingly more precise and interesting results from lattice QCD calculation of these 
rare decays in the future. Sufficiently accurate theoretical predictions of the form fac- 
tors are essential for making maximal use of experimental measurements, and may 
yield interesting constraints on new physics in the future. 
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Figure 2: The left panel is from C. Bouchard in the HPQCD collaboration. It is 
the preliminary result of i? — )■ Kl^l~ /+,o,t measured on the lattices. Only sta- 
tistical error is included and it is about 1% accuracy. The right panel is quoted 
from Ref. [211] • It is the preliminary result of the T12 in 5 — )■ K*l^l~ process from 
Cambridge/W&M/Edinburgh collaborations. 

6 Acknowledgements 

The author thanks the usefully discussions with the members in the FNAL/MILC col- 
laborations. The author also want to thank C. Bouchard and M. Wingate provide the 
progress of the works from the HPQCD collaboration and Cambridge/W&M/Edinburgh 
collaborations. 



References 

[1] B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 102, 091803 (2009) 
|arXiv:0807.4Tl9l [hep-ex] ] . 

[2] J. -T. Wei et al. [BELLE Collaboration], Phys. Rev. Lett. 103, 171801 (2009) 
[arXiv:0904.0770l [hep-ex] ] . 

[3] T. Aaltonen et al. [CDF Collaboration], Phys. Rev. Lett. 106, 161801 (2011) 
[arXiv:1101. 10281 [hep-ex] ] . 

[4] RAaij et al. [LHCb Collaboration], Phys. Rev. Lett. 108, 181806 (2012) 
[arXiv:1112.35T5l [hep-ex]]. 



8 



[5] RAaij et al. [LHCb Collaboration], Phys. Rev. D 85, 112013 (2012) 
[arXiv: 1202.62671 [hep-ex]]. 

[6] RAaij et al. [LHCb Collaboration], JHEP 1207, 133 (2012) [ arXiv:1205.3"422l 
[hep-ex]]. 

[7] J. L. Hewett, H. Weerts, R. Brock, J. N. Butler, B. C. K. Casey, J. Collar, A. de 
Govea and R. Essig et al., arXiv: 1205.2671 [hep-ex]. 



D. Becirevic, V. Lubicz and F. Mescia, Nucl. Phys. B 769, 31 
[hep-ph/0611295j . 

[9] L. Del Debbio et al. [UKQCD Collaboration], Phys. Lett. B 416, 392 
[hep-lat/9708008j . 

[10] A. Abad a et al. [APE Collaboration], Phys. Lett. B 365, 275 
[hep-lat /9503020|. 

[11] T. Bhattacharya and R. Gupta, Nucl. Phys. Proc. Suppl. 42, 935 
| hep4at/950T01 6|. 

[12] K. C. Bowler et al. [UKQCD Collaboration], Phys. Rev. Lett. 72, 1398 
|hep-lat/9311004| . 

[13] C. W. Bernard, P. Hsieh and A. Soni, Phys. Rev. Lett. 72, 1402 
[hep-lat/9311010] . 

[14] A. Al-Haydari et al. [QCDSF Collaboration], Eur. Phys. J. A 43, 107 
[arXiv:0903.1664 [hep-lat]]. 



[15] C. W. Bernard, T. Burch, K. Orginos, D. Toussaint, T. A. DeGrand, C. E. De- 



tar, S. Datta and S. A. Gottheb et al, Phys. Rev. D 64, 054506 
[hep-lat/0104002] . 



2007) 
1998) 
1996) 
1995) 
1994) 
1994) 
2010) 



2001) 



[16] R. Zhou, S. Gottlieb, J. A. Bailey, D. Du, A. X. El-Khadra, R. D. Jain, A. S. Kro- 
nfeld and R. S. Van de Water et a/.. arXiv:1211.1390 [hep-lat]. 

[17] R. Zhou et al. [Fermilab Lattice and MILC Collaborations], PoS LATTICE 2011, 
298 (2011) [arXiv:1111.0981l [hep-lat]]. 

[18] C. M. Bouchard, G. P. Lepage, C. J. Monahan, H. Na and J. Shigemitsu, 
Iar2£ivil2l0.6992 [hep-lat]. 

[19] Z. Liu, S. Meinel, A. Hart, R. R. Horgan, E. H. MuUer and M. Wingate, PoS 
LAT 2009, 242 (2009) |arXiv:0911. 23701 [hep-lat]]. 



9 



[20] Z. Liu, S. Meinel, A. Hart, R. R. Horgan, E. H. MuUer and M. Wingate, 
larXiv:1101.2726l [hep-ph]. 

[21] A. Bazavov, D. Toussaint, C. Bernard, J. Laiho, C. DeTar, L. Levkova, M. B. Ok- 
tay and S. Gottlieb et al, Rev. Mod. Phys. 82, 1349 (2010) jarXiv:090 3.3598. 
[hep-lat]]. 

[22] E. Follana et al. [HPQCD and UKQCD Collaborations], Phys. Rev. D 75, 054502 
(2007) |liep-lat/0610092|. 

[23] B. Sheikholeslami and R. Wohlert, Nucl. Phys. B 259, 572 (1985). 

[24] A. X. El-Khad ra, A. S. Kronfeld and P. B. Mackenzie, Phys. Rev. D 55, 3933 
(1997) I hep-lat /9604004| . 

[25] A. S. Kronfeld, Phys. Rev. D 62, 014505 (2000) |hep-lat/0002008;|. 

[26] M. B. Oktay and A. S. Kronfeld, Phys. Rev. D 78, 014504 (2008) 
[arXiv:0803.0523l [hep-lat]]. 

[27] G. P. Lepage, L. Magnea, C. Nakhleh, U. Magnea and K. Hornbostel, Phys. Rev. 
D 46, 4052 (1992) |he p-lat/9205007] . 

[28] R. R. Horgan, L. Khomskii, S. Meinel, M. Wingate, K. M. Foley, G. P. Lep- 
age, G. M. von Hippel and A. Hart et al, Phys. Rev. D 80, 074505 (2009) 
[jiX:iy:0906.0945 [hep-lat]] . 

[29] C. Bernard, M. Golterman and Y. Shamir, Phys. Rev. D 77, 074505 (2008) 
[arXiv:0712.256 [hep-lat]]. 

[30] M. Beneke, T. Feldmann and D. Seidel, Nucl. Phys. B 612, 25 (2001) 
|hep-ph/0106067| . 

[31] M. Beneke, T . .Feldmann and D. Seidel, Eur. Phys. J. C 41, 173 (2005) 
|hep-ph/ 0412400]. 

[32] K. G. Chetyrkin, M. Misiak and M. Munz, Phys. Lett. B 400, 206 (1997) 
[Erratum-ibid. B 425, 414 (1998)] |hep-ph/9612313j . 

[33] C. Bobeth, G. Hiller and D. van Dyk, JHEP 1007, 098 (2010) [ arXiv:1006.5013l 
[hep-ph]]. 

[34] J. A. Bailey, C. Bernard, C. E. DeTar, M. Di Pierro, A. X. El-Khadra, 
R. T. Evans, E. D. Freeland and E. Gamiz et a/., Phys. Rev. D 79, 054507 
(2009) [arXiv:0811.3640l [hep-lat]]. 



10 



[35] J. Foley, K. Jimmy Juge, A. O'Cais, M. Peardon, S. M. Ryan and J. -I. SkuUerud, 
Comput. Phys. Commun. 172, 145 (2005) |hep-lat/0505023). 

[36] G. P. Lepage, B. Clark, C. T. H. Davies, K. Hornbostel, P. B. Mackenzie, 
C. Morningstar and H. Trottier, Nucl. Phys. Proc. Suppl. 106, 12 (2002) 
[hep-lat /01101 75|. 

[37] K. Hornbostel, G. P. Lepage, C. T. H. Davies, R. J. Dowdall, H. Na and 
J. Shigemitsu, Phys. Rev. D 85, 031504 (2012) |arXiv:1111 .1363 [hep-lat]]. 

[38] C. Aubin and C. Bernard, Phys. Rev. D 76, 014002 (2007) [ arXiv: 0704. 0795) 
[hep-lat]]. 

[39] A. Bazavov, et al. [The MILG Collaboration], larXivil212J768l [hep-lat]. 

[40] S. Durr, Z. Fodor, C. Hoelbling, S. D. Katz, S. Krieg, T. Kurth, L. Lellouch and 
T. Lippert et al, Phys. Lett. B 701, 265 (2011) jarXivilOll. 24031 [hep-lat]]. 

[41] S. Durr, Z. Fodor, C. Hoelbling, S. D. Katz, S. Krieg, T. Kurth, L. Lellouch and 
T. Lippert et al, JHEP 1108, 148 (2011) |arXiv:1011.2"7ni [hep-lat]]. 

[42] J. A. Bailey et al [Fermilab and Lattice and MILG Collaborations], PoS LAT- 
TICE 2011, 270 (2011) |arXiv:1111.5"47n [hep-lat]]. 

[43] H. Na, C. T. H. Davies, E. Follana, G. P. Lepage and J. Shigemitsu, Phys. Rev. 
D 82, 114506 (2010) [arXiv: 1008.4562 [hep-lat]]. 



11 



Fpara NLO SU(2) ChPT fit. p-val=0.41 




0.02/0.05 
0.01/0.05 
0.007/0.05 
0.005/0.05 
0.0124/0.031 
0.0062/0.031 
0.0047/0.031 
0.0031/0.031 
0.00155/0.031 
0.0036/0.0188 
0.0018/0.0188 
cont 



1 ^ ' ' ' ' ' 1 

0.8 1 1.2 1.4 1.6 1.8 2 

(n Ek) 



Fperp NLO SU(2) ChPT fit. p-val=0.41 




0.02/0.05 
0.01/0.05 
0.007/0.05 
0.005/0.05 
0.0124/0.031 
0.0062/0.031 
0.0047/0.031 
0.0031/0.031 
0.00155/0.031 
0.0036/0.0188 
0.0018/0.0188 
cont 



(n Ek) 



